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O ABSTRACT 

o 

' A study of the ultra high energy neutrino detection performances of a km 3 Neu- 

CL|! trino Telescope sitting at the three proposed sites for "ANTARES" , "NEMO" 

and "NESTOR" in the Mediterranean sea is here performed. The detected 
charged leptons energy spectra, entangled with their arrival directions, provide 
an unique tool to both determine the neutrino flux and the neutrino-nuclcon 
cross section. 



Neutrinos are one of the main components of the cosmic radiation in the ultra- 
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high energy (UHE) regime. Although their fluxes are uncertain and depend on the 
production mechanism, their detection can provide information on the sources and 
origin of the UHE cosmic rays. 

From the experimental point of view the detection perspectives are stimulated by 
the several proposals and R&D projects for Neutrino Telescopes (NT's) in the deep 
water of the Mediterranean sea, namely ANTARES NESTOR 2) an d NEMO^, which in 
the fut ure c ould lead to the construction of a km 3 telescope as pursued by the KM3NeT 
project Actually, on the ANTARES site, a smaller telescope with a surface area 
of 0.1 km 2 is already under constructi on ® . A further project is IceCube, a cubic- 
kilometer under-ice neutrino detector EBE^ currently being deployed in a location 
near the geographic South Pole in Antarctica. IceCube applies and improves the 
successful technique of AMANDA to a larger volume. 

Although NT's were originall y thought as detector s, their capability as v T 
detectors has become a hot topic H0 | 1 1 | 12 | 13 | 14 | 15|16 | 1 7 | 18t in view of the fact that 

flavor neutrino oscillations lead to nearly equal astrophysical fluxes for the three 
neutrino flavors. Despite the different behavior of the produced tau leptons with 
respect to muons in terms of energy loss and decay length, both and v T detection 
are sensitive to the matter distribution near the NT site. Thus, a computation of the 
event detection rate of a km 3 telescope requires a careful analysis of the surroundings 
of the proposed site. The importance of the elevation profile of the Earth surface 
around the detector was already found of some relevance in Ref. where s ome of 
the present authors calculated the aperture of the Pierre Auger Observatory I 20 | 21 | 
for Earth-skimming UHE z/ r 's. Indeed, air shower experiments can be used as NT's 
at energies > 10 18 eV, a topic recently reviewed in ^2}. In Ref. ^2} it is estimated 
the effective aperture for v T and detection of a km 3 NT in the Mediterranean 
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Figure 1: The surface profile of the area near the ANTARES site (red spot) at 42° 30 N, 07° 00 E. 
The black curve represents the coast line. The sea plateau depth in the simulation is assumed to 
be 2685 m. The effective volume starts at an height of 100 m from the seabed, to account for the 
spacing of the first photomultipliers as foresee by the current designs. 




Figure 2: The surface profile of the area near the NEMO site (red spot) at 36° 21 N, 16° 10 E. The 
black curve represents the coast line. The sea plateau depth used in the simulation is 3424 m. The 
effective volume starts at an height of 100 m from the seabed, to account for the spacing of the first 
photomultipliers as foresee by the current designs. 




Figure 3: The surface profile of the area near the NESTOR site (red spot) at 36° 21 N, 21° 21 E. The 
black curve represents the coast line. The sea plateau depth in the simulation is assumed to be 4166 
m. The effective volume starts at an height of 100 m from the seabed, to account for the spacing of 
the first photomultipliers as foresee by the current designs. 



sea placed at any of the three locations proposed by the ANTARES, NEMO and NESTOR 
collaborations. The characteristics of the three site surface profiles are compared 
by using the DEM of the different areas. 

In the present paper we further develop the approach of Ref. ^3 j n order to apply 
the detection of UHE v as a tool to simultaneously measure the UHE neutrino flux 
and the z/-N cross section in extreme kinematical regions. 

Following the formalism developed in we define the km 3 NT fiducial volume 
as that bounded by the six lateral surfaces S a (the subindex a=D, U, S, N, W, and 
E labels each surface through its orientation: Down, Up, South, North, West, and 
East), and indicate with Q a = (9 a ,4> a ) the generic direction of a track entering the 
surface E . The scheme of the NT fiducial volume and two examples of incoming 
tracks are shown in Fig. HI We introduce all relevant quantities with reference to v T 
events, the case of being completely analogous. 

Let d<f>„/(dE u dQ a ) be the differential flux of UHE v r + v T . The number per unit 
time of r leptons emerging from the Earth surface and entering the NT through S a 
with energy E T is given by 

f^r) = S dQa S dSa I dK d dE^ a) / dEr cos {6a) K{Ev ' Et] fa> Qa) • (1) 




Figure 4: The angle definition and the fiducial volume of a km 3 NT. 

The kernel k T a (E v , E T ; f a , Q a ) is the probability that an incoming u T crossing the 
Earth, with energy E v and direction Q a , produces a r-lepton which enters the NT 
fiducial volume through the lateral surface dS a at the position r a with energy E T (see 
Fig.HJfor the angle definition). If we split the possible events between those with track 
intersecting the rock and the ones only crossing water, the kernel k T a {E v , E T ; f a , Q a ) 
is given by the sum of these two mutually exclusive contributions, 

k T a (E v , E T ■ f a , tt a ) = k T /(E u , E T ; r a , o ) + k T a ' w (E v , E T ; f a , fi a ) . (2) 

For an isotropic flux we can rewrite Eq. (jTJ), summing over all the surfaces, as 

dNl r ' w = ^[dE u [ dE T I dQ a I dS a fi- d ^ E K) ) cos (6 a ) kl^\E v , E T ; r , fi a ) 



By using this expression one can also define the total aperture A T( - r > w ^(E u ), with "r" 
and 'W denoting the roc/c and water kind of events, respectively, 



dt 



dE„ 



' 1 d$ v (E v ) 
4vr dE y 



(4) 



where 



A T < r ' tt >(£7„) = E / d ^ / dfi a | dS a cos (0 a ) fc^^ , £ r ; r a , fi a ) . (5) 



(3) 



Of course, the same quantities can be defined for muons coming from the charged- 



current interactions of u^. 



In Fig. [5] we compare the detection performances of a km 3 NT placed at one of the 
three sites in the Mediterranean sea. The NESTOR site shows the highest values of the 
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Figure 5: A comparison of the effective apertures A T ^ r,w ^\E v ) for the three NT sites. We 
plot the ratios [A T « r ' w » (NESTOR) - A T (( r > tu ')(NEM0)]/ J 4 r (( r > to ')(NEM0) and [A T « r ^» (ANTARES) - 
A T ^ r ' w ^(mK0)]/A T ^ r ' w ^(mm) versus the neutrino energy. 



r-aperture for both rock and water, due to its larger depth and the particular matter 
distribution of the surrounding area, while the lowest rates are obtained for ANTARES. 
The aperture in the three sites can be quite different at high energy but, in order to 
get the expected number of UHE events per year, one has to convolve the aperture 
with a neutrino flux which typically drops rapidly with the energy. Although the 
percentage value of the matter effects remains unchanged, in this very low statistics 
regime they can be hardly distinguished; still, they can be enhanced by an appropriate 
choice of the detector shape. 

Knowing the aperture of the NT at each site, we can compute the expected r event 
rate, once a neutrino flux is specified. In Table [1] these rates are shown assuming a 



Surf. 


ANTARES 


NEMO 


NESTOR 


D 


0.0059/0 


0.0059/0 


0.0058/0 


U 


0/0.1677 


0.0002/0.2133 


0.0002/0.2543 


s 


0.0185/0.1602 


0.0256/0.1773 


0.0240/0.2011 


N 


0.0241/0.1540 


0.0229/0.1823 


0.0321/0.1924 


W 


0.0212/0.1584 


0.0335/0.1691 


0.0265/0.2002 


E 


0.0206/0.1589 


0.0190/0.1875 


0.0348/0.1907 


Total 


0.090/0.799 


0.107/0.929 


0.123/1.039 



Table 1: Estimated rate per year of rock/water t events at the three km 3 NT sites for a GZK-WB 
flux 23 '' 25 '. The contribution of each detector surface to the total number of events is also reported. 
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Figure 6: In the two panels are reported the number of events (/i + r) collected in three years from 
a km 3 NT. The left panel concerns the events with energy lost in the detector in the range 10 5 -10 8 
GeV, whereas on the right the events have an energy deposited larger than 10 s GeV. See the text 
for further details. 

GZK-WB flux 1^123}. The effect due to the local matter distribution is responsible 
for the N-S, W-E and NE-SW asymmetries for the ANTARES, NEMO and NESTOR sites, 
respectively, as expected from the matter profiles shown in Figs. [TJ El and O These 
matter effects, for the specific UHE flux considered (GZK-WB), correspond to an 
enhancement of rock events which goes from 20 to 50% for the three sites, respectively, 
and a screening factor for water events from 3 to 10%. The largest relative difference 
among lateral surfaces is in the case of W/E for NEMO, where the huge wall to the west 
of the site (see Fig. [2]) improves the rate by about 75%, almost a factor 2! Notice 
also that the water events from the U surface are basically proportional to the depth. 

Due to the dependence of Eq. ([3]) on the neutrino flux and the different behavior 
of kl' r (E u , E T ; f a , Q a ) and k T ^ w {E v , E T ; r a , Q a ) as functions of the neutrino-nucleon 
cross section, o~qq, one can imagine to use the detected events, properly binned for 
energy loss and arrival direction, in order to obtain information on both the neutrino 
flux and the neutrino-nucleon cross section. In particular, since the real observable 
is the energy deposited in the detector and not the energy and/or the nature of the 
charged lepton, either /i or r, crossing the NT, one must sum the two contributions. 
In fact, the events whose topology allows for determining the nature of the charged 
lepton are a negligible fraction of the expected total number. 

In the two panels of Figure [6] are reported the number of events (/x + t) collected in 
three years from a km 3 NT. In particular the left panel concerns the events where the 
energy deposited in the detector is in the range 10 5 -10 8 GeV, whereas the one on the 
right reports events where the energy lost is larger than 10 8 GeV. For each panel the 
events have been split in three bins according to their arrival directions (0° represents 
the vertical downgoing direction). Fixing the panel and the arrival direction range, 
the three bars of the histogram represents three different neutrino fluxes and o v q G 
chosen. In particular from left to right we have the GZK-WB 23 )' 24 ) flux and standard 
cross section, the GZK-WB 23 )' 24 ) flux and three times the standard cross section and 
finally, the more copious GZK-H 23 )' 24 ) flux and standard cross section. Note that 
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Figure 7: The number density of yearly events (/z plus r) as function of the energy deposited by 
the charged lepton AEi for different values of A and with B = and C = D = 1 . The solid line 
represent the standard scenario. 
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Figure 8: The number density of yearly events (/z plus r) as function of the energy deposited by 
the charged lepton AEi for different values of B and with A = 3 and C = D = 1. The solid line 
represent the standard scenario. 
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Figure 9: The number density of yearly events (fj. plus r) as function of the energy deposited by 
the charged lepton AEi for different values of D and with A = C = 1 and B = 0. The solid line 
represent the standard scenario. 



the total number of events for GZK-WB with 3 o v q C and for GZK-H with standard 
cross section are the same. The different flux/cross section configurations can be 
disentangled by observing the different behavior of the height of the bars as function 
of the energy lost and arrival direction. 

In order to study the sensitivity to both neutrino flux and Oq C it is necessary to 
parameterize their standard expression and the possible departure from it. In partic- 
ular by using a standard Waxman-Bahcall ^ (C=D=1) as a conservative reference 
for the neutrino flux we allow for a variation of its steepness via the exponent D and 
for the normalization through the multiplicative factor C, namely: 

(f) WB = 1.3 ■ 1(T 8 C ef D GeV" 1 cm" 2 s" 1 sr" 1 . (6) 

In the same way, for the neutrino-nucleon cross section one can parameterize the 
presence of new physics by assuming a departure from the standard expression j n 
terms of two free parameters, A and B, whose standard values are A = 1 and B = 0: 

' 0.677 e°- 492 ; 2.00 ■ 10 4 < e v < 1.20 • 10 7 

1.20 ■ 10 7 <e v < 1.20 • W +B (7) 

5.54 ■ l 0.363(l-A)(7.08+B) e 0.363A . 1.20 ■ 10 7+B < € v 

where e v = E v j GeV. In particular B fixes the energy value where new physics appears 
and A is the change in the energy slope of cr^c- m Figures [3, [8] and [9] it is reported 
the effect of the variation of the single parameter A, B, D on the number density of 
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yearly events (/i plus r) as a function of the energy deposited by the charged lepton 
AEi. The factor C has been fixed to its standard value (C = 1) since it is just a 
normalization and thus simply correlated to the exposure time needed to achieve the 
proper event statistics. 

The quite relevant effect shown by Figures [7J [8] and [9] supports once more the 
idea that a km 3 NT can provide a real chance to both measure UHE neutrino flux 
and the neutrino-nucleon cross section in the extreme kinematical region, where some 
new physics could appear. Of course the real feasibility of such measurements will 
crucially depend of the size of the neutrino flux which fixes the time required to 
reach a reasonable statistics. Then, while in this exercise we adopted an extreme 
conservative point of view working with Waxman-Bahcall like fluxes, one can wish 
for a more optimistic real situation. 
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